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FAR-ULTRAVIOLET & Ha SPECTROSCOPY OF SNR 0057-7226 IN THE SMC H II REGION N66 
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ABSTRACT 

N66 /NGC 346 is the largest and brightest H II region in the Small Magellanic Cloud and contains at 
least one known supernova remnant SNR 0057— 7226. Optical emission from the remnant is overwhelmed 
by the bright photoionized emission from the nebula, but the remnant has been detected by way of far 
ultraviolet absorption lines. Here we present data from the Far Ultraviolet Spectroscopic Explorer 
(FUSE) satellite showing strong O VI and C III emission from a position at the edge of SNR 0057— 7226. 
We also present high-resolution, long-slit Ha spectra across N66 showing high- and low- velocity emission 
corresponding closely to the X-ray boundaries of the supernova remnant. We use these FUV and optical 
data to determine the physical parameters of the shock and interaction geometry with N66. We find 
that ionizing photons from the many massive cluster stars nearby likely affect the ionization balance 
in the post-shock gas, hindering the production of lower-ionization and neutral species. We discuss the 
importance and inherent difficulty of searching for supernova remnants in or near bright H II regions 
and suggest that the far ultraviolet provides a viable means to discover and study such remnants. 

Subject headings: H II regions— ISM:individual(N66, SNR 0057-7226, NGC 346)— Magellanic 
Clouds — supernova remnants 



1. INTRODUCTION 

The giant H II region N66 is the largest, brightest star 
forming region in the Small Magellanic Cloud (SMC) . The 
populous star cluster NGC 346 embedded in the nebula is 
home to several dozen early type stars including a dozen 
of spectral types 07 and earlier. Two arcminutes to the 
east lies the massive, luminous WN+OB system HD 5980 
which underwent a luminous blue variable (LBV) outburst 
in 1994 (Koenigsberger et al. 1995; Barba et al. 1995). 
These stars are the main source of ionization in the neb- 
ula. 

N66 also contains at least one known supernova rem- 
nant (SNR), SNR 0057-7226. Einstein X-ray observations 
showed a bright source in N66 tentatively identified as a 
SNR (Inoue, Koyama, & Tanaka 1983; Wang & Wu 1992). 
Radio observations by Mills et al. (1982) showed bright, 
extended, thermal emission at 408, 843 and 5000 MHz 
corresponding to the H II region, but no clear sign of the 
SNR. Optical observations also showed only the extremely 
bright nebular emission but no filamentary structure asso- 
ciated with the SNR. However, Chu & Kennicutt (1988) 
discovered faint, high- velocity Ha emission in a long-slit 
spectrum across N66 and suggested the existence of a SNR. 
Subsequently, Ye, Turtle, & Kennicutt (1991) subtracted 
smoothed Ha data from 843 MHz radio continuum obser- 
vations and revealed a remnant —3.2' (~55 pc at an SMC 
distance of 59 kpc) in diameter and confirmed the remnant 
SNR 0057-7226. X-ray spectroscopy of SNR 0057-7226 
with ASCA (Yokogawa et al. 2000) showed that the spec- 
trum can be modeled as a hot thermal plasma. 



High-resolution X-ray observations of N66 with the 
Chandra X-ray Observatory (Naze et al. 2002) show a 
130" x 100" region of extended, center-bright, thermal 
emission corresponding to the SNR. These high-resolution 
data show SNR 0057— 7226 to have relatively uniform sur- 
face brightness and no clear temperature gradient from 
center to rim. Difference imaging techniques applied to re- 
cent radio observations (Filipovic et al. 2003) at 1.42, 2.37, 
and 4.80 GHz (20, 13 and 6 cm) show SNR 0057-7226 to 
be a non-thermal, limb-brightened shell with a spectral in- 
dex of approximately a—~ 0.17 (S„ oc v a ). The X-ray and 
radio emitting regions are similar in extent at a resolu- 
tion of -20". Hence, the Ye, Turtle, & Kennicutt (1991) 
estimate of angular size is probably too large, and the 
X-ray/radio SNR is roughly 36 x 28 pc in extent. Thus, 
SNR 0057— 7226 may be a member of the "mixed morphol- 
ogy" class of SNRs (Rho & Petre 1998), with a radio shell 
and filled-center X-ray emission, although the presence of 
HD 5980 along the sight line confuses the situation. Naze 
et al. (2002) resolve X-ray emission from HD 5980 itself for 
the first time, and note similarities with the galactic object 
rj Carinae, which has diffuse surrounding X-ray emission 
(Seward & Mitchell 1981). 

Naze et al. (2002) mention a possible association of 
HD 5980 and SNR 0057-7226, but FUV absorption stud- 
ies have clearly determined that HD 5980 is behind the 
SNR. FUV observations toward HD 5980 show the ex- 
pected Galactic and SMC absorption at vi sr w0 km s _1 
and vi sr «+150 km s _1 , respectively, as well as high- 
velocity absorption indicating the receding shell of an in- 
tervening SNR. Absorption at +300 km s _1 was found in 
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IUE spectra by de Boer & Savage (1980); this was later 
confirmed by Fitzpatrick & Savage (1983) who suggested 
an intervening SNR. Koenigsberger et al. (2001) detected 
absorption systems at +300 and +330 km s _1 in C IV, 
N V and several other ions in high-resolution HST STIS 
spectra. They also see very weak components at +20 
and +50 km s -1 , possibly consistent with the approach- 
ing side of the remnant if indeed these features are real. 
Hoopes et al. (2001) also see absorption at ~+300 km s _1 
in C III and O VI in FUSE observations of HD 5980 and 
conclude that the star must lie behind the receding shell 
of SNR 0057-7226. Hoopes et al. (2001) also find that a 
FUSE spectrum of Sk80 (AV232), an 071a star -1' south- 
east of HD 5980 but also seen in projection within the X- 
ray emission, shows none of the high-velocity absorption 
seen toward HD5980. Since Sk80 lies closer to the X-ray 
edge of the remnant, we might expect a lower line-of-sight 
velocity for SNR gas in this sight line. However, the O VI 
absorption profile near the SMC systemic velocity (^150 
km s _1 ) is indistinguishable from that of HD 5980. Thus, 
Sk80 must either lie within or in front of the remnant. 

Far-ultraviolet (FUV) observations are useful in the 
study of SNRs. Several strong emission lines including 
C IV AA1548,1550, N VAA1238,1242, and especially O VI 
AA1032,1038 are good shock diagnostics and are produced 
in regions where the temperature is lower than typical X- 
ray producing regions and higher than regions of bright op- 
tical emission (Sutherland & Dopita 1993). O VI emission, 
in particular, arises in gas at temperatures near 3 x 10 5 K- 
a condition almost never reached by photoionization-and 
is thus a sensitive tracer of shock-heated material. 

In this paper we present a FUSE observation of emis- 
sion at the X-ray edge of SNR 0057-7226 that shows 
strong emission in O VI AA1032,1038 and C III A977. 
We also present new longslit echelle observations showing 
high- and low-velocity material in Ha across the face of 
SNR 0057— 7226. The observations are described in §2. In 
§3, we discuss the SNR kinematics as revealed by the opti- 
cal and FUV data. In §4, we compare the O VI, C III, and 
Ha measurements with the predictions of shock models to 
derive shock velocity, preshock density and ram pressure at 
one location in the SNR. Using the same models with the 
FUV absorption data of Hoopes et al. (2001) and Koenigs- 
berger et al. (2001), we derive similar values for material 
on the sight line toward HD 5980. We propose a physi- 
cal picture for SNR 0057-7226 in relation to N66. Finally 
we discuss the implications these observations have for the 
detection of SNRs in bright nebulae and OB associations. 
A summary is presented in §5. 

2. OBSERVATIONS 

Figure 1 shows images of SNR 0057— 7226 and its sur- 
roundings. The left panel shows an Ha+[N II] emission- 
line image taken with the Curtis Schmidt Telescope at 
CTIO in 1999, December. These data are described more 
thoroughly by Danforth et al. (2002). The right panel 
shows a recent X-ray image of NGC 346 from the Chandra 
X-ray Observatory, obtained courtesy of Naze et al. (2002). 
These data show SNR 0057— 7226 clearly as a bright, ex- 
tended source. Unfortunately, the ACIS detector gaps par- 
tially obscure the southwestern edge of the SNR and some 
of the core of NGC 346 itself. 



2.1. FUSE Observations 

Six stars within N66 /NGC 346 are part of the FUSE Hot 
ISM Team program and are marked with small boxes in 
Figure 1. Four of these targets are within the crowded 
core of the cluster itself. The cluster targets were ob- 
served with the 4"x20" MDRS aperture. The remaining 
two (HD 5980 and Sk80) lie about two arcminutes to the 
east of the cluster in projection toward SNR 0057— 7226 
and were observed with the 30"x30" LWRS aperture. 

FUSE has three spectral apertures (Moos et al. 2000); 
the LWRS and the MDRS apertures are centered 209" 
from each other. Thus, while the MDRS aperture was 
being used to observe stars in NGC 346, the LWRS aper- 
ture fell serendipitously on other parts of the sky. In one 
case (P20305, henceforth Position 1), the LWRS aper- 
ture fell 75" east of HD 5980 on the eastern X-ray edge 
of SNR 0057— 7226. This spectrum shows bright emission 
lines of C III A977.020 and O VI AA1031.926, 1036.617. 
A very faint continuum underlies the bright emission lines 
at a flux level of ~ 10~ 14 erg s _1 cm~ 2 A -1 . In another 
case (P20302, henceforth Position 2) the LWRS aperture 
fell in a region north of HD 5980 outside the boundaries of 
the SNR's X-ray emission. This sight line shows a weak 
continuum spectrum similar to that at Position 1. We use 
this sight line to characterize the FUV background of N66. 
The FUSE observations are summarized in Table 1. Fig- 
ure 1 shows the positions of the FUSE LWRS apertures 
for Positions 1 and 2 as large boxes. 

The third FUSE aperture, the l"25x20" HIRS aper- 
ture, lies halfway between the other apertures. In several 
observations of stars in NGC 346 the HIRS aperture fell on 
X-ray bright regions of SNR 0057-7226. However, no flux 
was detected in these spectra. Since the HIRS aperture 
has less than 3% of the collecting area of LWRS for de- 
tecting diffuse emission, this negative result is consistent. 

FUSE has four independent spectral channels (in FUSE 
nomenclature, LiFl, LiF2, SiCl, and SiC2), and because 
of thermal effects onboard the spacecraft, relative chan- 
nel alignments shift in a complicated manner (Moos et al. 
2000; Sahnow et al. 2000). This can cause significant 
complications for diffuse source observations (cf. Blair, 
Sankrit, & Tulin 2002). The spacecraft guides on the LiFl 
channel, which also has the highest FUV throughput near 
1035A. The other three channels tend to drift as much as 
8" over the course of a few orbits with respect to LiFl if 
no realignments are made. For this reason, we use only 
LiFl data for stellar targets whenever possible. LiFl does 
not cover wavelengths below 987A. Thus we use the lower- 
throughput SiC2 channel to observe C III A977. All four 
channels cover the 1035A region containing O VI, so a di- 
rect comparison of C III and O VI is possible. The six 
stellar datasets were reduced using the methods described 
by Danforth et al. (2002) and references therein. Spectral 
resolution in these data is 15-20 km s _1 . 

The LWRS data were reduced using CALFUSE v2.0.5 
and observations were combined by channel. The region 
around 1032A was examined in all four channels to deter- 
mine the relative shifts that exist in the wavelength so- 
lutions relative to LiFl. We increased the signal-to- noise 
ratio by coadding the two LiF channels to create a com- 
posite O VI spectrum. This does not significantly degrade 



3 



the spectral resolution^ which for an aperture-filling ex- 
tended source is ~0.34A (106 km s _1 at 1032A). O VI flux 
is present in the two SiC channels as well and closely re- 
sembles that in the LiF channels, but with lower S/N. The 
two SiC channels were coadded for the C III emission. Fig- 
ure 2 shows these data in comparison to Position 2 which 
sampled a sight line outside the SNR. 

2.2. CTIO 4-m Echelle Spectra 

Seven longslit echelle spectra of the N66 region were 
obtained on the 4-m Blanco telescope at the Cerro Tololo 
Interamerican Observatory (CTIO) during 2000 December 
and 2001 December. These observations are summarized 
in Table 2 and the slit positions are shown in Figure 1 as 
dashed lines. For both observing runs, a 79 line mm -1 
echelle grating and long-focus red camera were used. The 
decker was opened to its maximum extent giving a useful 
slit length on the sky of 3f5. The slit width was set to 1'.'65. 
The cross-dispersing grating was replaced with a flat mir- 
ror and an order-separating filter (A C =6580A, AA=150A) 
was used to isolate the spectral region around the Ha line. 
[N II] AA6548,6583 also falls within the spectral range sam- 
pled. 

Data reduction was performed with IRAF 5 using the 
QUADPROC and longslit packages. Trim, bias and dark 
current calibration were performed first. The images were 
next processed to remove the tilt of the CCD axis with 
respect to the slit. This rendered the stellar spectra hor- 
izontal along a row. Thorium- Argon lamp spectral lines 
were traced and fit to correct for curvature in the cross- 
dispersion direction and provide a dispersion solution. Fi- 
nal wavelength calibration was performed using the terres- 
trial airglow lines present in the observations themselves. 
The result is a rectified spectrum with dispersion along 
the x-axis and the spatial dimension along the y-axis. The 
calibrated data cover ~3.'5 of sky at a spatial resolution of 
<~1" (0'.'267 pix -1 ) and spectral resolution of ~10 km s _1 
(0.082A pix" 1 = 3.7 km s" 1 pix -1 ). We present three of 
the spectra in two-dimensional format in Figure 3, where 
the wavelength scale has been converted to velocity. 

3. ANALYSIS 

3.1. FUV Emission 

The LWRS aperture of Position 1 lies 75" east of 
HD 5980 and coincides with the edge of bright, soft X-ray 
emission (Figure 1). Strong O VI emission is seen in both 
lines of the doublet as well as C III A977.02 (Figure 2). 
We detect no other lines in the FUSE wavelength range 
(905-1 187A) besides airglow. In addition to the emission 
lines, a faint continuum is seen in the Position 1 spec- 
trum (Fa ~ 10~ 14 erg s _1 cm -2 A -1 ). The spectrum is 
consistent with an early-type star and shows many of the 
same ISM absorption features seen in other FUSE spec- 
tra. The unidentified star in the LWRS aperture visible in 
the optical image (Figure 1) is not the source of the FUV 
continuum. It does not appear in FUV UIT (A = 1550A) 
images of the region (Parker et al. 1998) and is of spectral 
type A or later. 

In contrast to Position 1, there are no O VI and C III 
emissions detected at Position 2. This aperture fell an 

5 IRAF is distributed by the National Optical Astronomy Observatories, operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 



arcminute northeast of the edge of the SNR in a field 
empty of bright stars, UV-bright or otherwise. How- 
ever, a similar stellar continuum can be seen in this spec- 
trum. The faint continuum observed at both positions is 
due to dust-scattered starlight from the many FUV-bright 
sources nearby. A similar effect at longer UV wavelengths 
was seen in IUE spectra of other regions in the SMC (Blair 
et al. 1989). 

Any emission from the SNR will suffer absorption from 
the intervening SMC and Galactic ISM. Fortunately, there 
is little molecular hydrogen absorption in the vicinity of 
1032A; the only ISM absorption comes from SMC and 
Galactic O VI. We use the continuum-normalized O VI 
absorption profile of Sk80 (Danforth et al. 2002, their Fig- 
ure 80) as our absorption model. A Gaussian profile di- 
vided by the absorption model was convolved with a 106 
km s _1 aperture function. The Gaussian parameters were 
varied and the resulting models were fit to the observed 
A1032 line profile. The best fit parameters were: vi sr =228 
km s-\ FWHM=159 km s"\ and peak flux=2.23x 10~ 13 
erg s _1 cm~ 2 A -1 . The total flux of the intrinsic Gaus- 
sian function is 1.27xl0 -13 erg s _1 cm~ 2 compared with 
an observed flux of 8. 68x10" 14 erg s _1 cm~ 2 . Both the fit 
to the data and the intrinsic Gaussian profile are shown in 
the top panel of Figure 4. 

The O VI A 1032 line should be twice the intensity of the 
A1038 line in the absence of resonance scattering and self- 
absorption. However, this ratio can vary as a function of 
velocity within the line profile. The observed A1032/A1038 
ratio is close to two on the blue-shifted side of the line but 
becomes significantly higher on the red-shifted wing. Since 
SMC and Galactic O VI absorption affect both lines in a 
proportional manner, we conclude that resonance scatter- 
ing and self-absorption is minimal but that significant ab- 
sorption from molecular hydrogen — 5-0(Pl)H 2 A1038.156 
and 5-0(R2)H 2 A1038. 690— affects the weaker line. 

We can characterize the A1038 emission using the in- 
trinsic emission profile determined from the A1032 line. 
In this case, the absorption model was allowed to vary; a 
simulated molecular absorption spectrum was calculated 
with the FUSE spectral simulator code fsim and the Sk80 
A1032 O VI absorption profile (scaled by a square root) 
was multiplied in. (Since T1038 = 5T1032 in the optically 
thin limit, 1(1038) = I Q exp(- ±71032) = I(1032) 1 /2.) The 
resulting model was multiplied by the intrinsic Gaussian 
profile determined above and convolved with the aperture 
function to simulate the observed A1038 line. 

We first fit the A1038 emission using two molecular hy- 
drogen components: a Galactic component (Shull et al. 
2000) (v=0 km s"\ N(H 2 )=1.6 x 10 16 cur 2 , b=5 km s"\ 
T=150 K) and an SMC component (v=115 km s" 1 , 
N(H 2 )=2 x 10 15 cm" 2 , b=15 km s"\ T=400 K) slightly 
modified from Tumlinson et al. (2002). However, this did 
not produce enough absorption on the red side of the emis- 
sion. A third component, representing swept-up material 
from the SNR at v~155 km s -1 improves the fit (Figure 4, 
middle panel), but requires an unphysically large turbu- 
lent velocity for a good fit. There could be many compo- 
nents at a range of velocities in the turbulent, post-shock 
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gas which might account for the enhanced absorption. No 
such high- velocity components are seen by Tumlinson et al. 
(2002). 

The C III A977 emission is strongly affected by reso- 
nance scattering and self-absorption. The FUSE spectrum 
shows that the Galactic and SMC components are satu- 
rated (Danforth ct al. 2002). Therefore we do not at- 
tempt to fit the observed C III profile. The measured flux, 
6.3xl0~ 14 erg s _1 cm~ 2 , is a lower limit to the intrinsic 
flux. Overlaying the derived O VI A1032 profile (Figure 4, 
bottom panel) provides a first order assessment of how this 
line might be affected if its intrinsic line width is similar 
to O VI. 

Apart from O VI and C III (and geocoronal airglow), 
there arc no other lines detected. In a spectrum of the 
Cygnus Loop, a typical middle-aged SNR, the next bright- 
est lines to O VI and C III were [Ne V] A1146 and S IV 
AA1063,1074 at -5% of the observed O VI flux (Blair, 
Sankrit, & Tulin 2002). Such weak emission would not 
be detected in the spectrum of SNR 0057— 7226, especially 
considering the lower abundances of these elements in the 
SMC. 

3.2. Ha and [N II] Emission 

The seven echelle observations cut across much of the 
N66 nebulosity in Figure 1 including the core of NGC 346 
and the outlying regions. Three spectra cross significant 
portions of the SNR and show faint high- and low-velocity 
emission. The data show both Ha and [N II] emission at 
a velocity resolution of ^15 km s _1 . 

We find the nebula itself to be surprisingly quiescent on 
the large scale given the violence of strong stellar winds 
and ongoing star formation. The bright nebular emis- 
sion can be fit with a single-component Gaussian with 
v lsr =157.7 km s _1 , FWHM=34.7 km s^ 1 . Assuming a 
temperature of 10 4 K in the photoionized gas, we derive a 
bulk turbulence of 23.8 km s _1 (b=14.3 km s" 1 ). This de- 
gree of turbulence is typical of Magellanic H II regions but 
much lower than that seen in giant H II regions such as 30 
Doradus and nebulae in M101 and M33 (Chu & Kennicutt 
1994; Smith & Weedman 1973, 1971, 1970). The thermal 
FWHM of [N II] is narrower than that of Ha by the square 
root of the ratio of their ionic masses (i.e. by \/l4=3.7). 
Thus [N II] emission is a more sensitive probe of veloc- 
ity structure than Ha. However, these lines are usually 
much weaker than Ha, particularly in the low-nitrogen 
SMC ISM, and are often completely undetectable in SMC 
echelle spectra. Four of the echelle spectra (both HD 5980 
pointings, NGC 346-WB3+4 and NGC 346- WB6) have us- 
able [N II] data. These show vi sr = 158.5 km s _1 and a 
FWHM=22.8 km s^ 1 . Assuming a temperature of 10 4 
K, we derive a bulk turbulence of 17.5 km s _1 (b=10.5 
km s _1 ), also typical for ordinary Magellanic H II regions. 

The two spectra which cut through NGC 346 show an 
expansion feature bounded by the bright Ha rim visible 
in Figure 1 surrounding the cluster core. (NGC346-WB1 
samples the same region as the latter two spectra, but has 
a high sky brightness masking any [N II] emission.) The 
data show a maximum expansion velocity of 14 km s _1 . 
This expansion is presumably driven by the stellar winds 
of the stars located within the bubble. The NGC 346 ex- 
pansion feature is also detected in the Ha data, but only 



in that a pair of Gaussians are needed to fit the broad 
Ha profile, rather than a single component. Meanwhile, 
the east- west HD 5980 data cut through N66C, a compact 
H II region to the west of HD 5980 and north east of the 
cluster core. The [N II] data reveal an expansion feature 
here with v cxp m 13 km s" 1 . 

High- and low-velocity Ha emission appear in the four 
echelle observations which cross SNR 0057- 7226 (HD 5980 
NS and EW, Sk80 NS and EW). Three of these spectra are 
shown in Figure 3; the fourth spectrum (Sk80 EW) is of 
low signal-to-noise and intersects only the edge of the rem- 
nant. The faint, higher velocity emission in these spectra 
arises in shocks associated with the SNR. It is weaker than 
the photoionized nebular emission by a factor of ~70 and 
has a very low signal to noise ratio. The SNR emission is 
quite patchy: blobs with widths of ~10" (~3 pc) are seen 
in the position- velocity diagrams (Figure 3). 

We find high-velocity features all across the face of 
SNR 0057—7226 with higher velocity emission typically to- 
ward the center suggesting an expanding shell. The fastest 
emission is at v~ +335 km s _1 about 15" (4 pc) south of 
HD 5980 (this is also the region of peak X-ray brightness). 
If the systemic velocity of the SNR is taken to be the same 
as the N66 ISM, this gives an expansion velocity of ~ 175 
km s _1 , consistent with strong O VI production by the 
shock front. 

The Ha emission from the approaching side of the SNR 
is also seen in Figure 3, though it is even fainter than the 
receding side. In the few cases where the blue-shifted ma- 
terial is discernible from the Gaussian wings of the bright 
nebular emission, we see velocities as low as vi sr ~ +50 
km s _1 . As with the receding side of the SNR, the mate- 
rial with the largest velocity offsets tends to be closest to 
the center of the remnant. 

Approximate spectrophotometry was performed on the 
echelle data using a flux-calibrated Ha+[N II] image of 
N66 obtained from Kennicutt & Hodge (1986). Comparing 
profiles from the image with wavelength-integrated inten- 
sity along the slit gave a conversion factor of ^4.5xl0~ 18 
erg s _1 cm~ 2 count -1 . We use this value to measure fluxes 
and surface brightnesses for different patches of SNR emis- 
sion in the spectra (Figure 3). A Gaussian profile was 
fit to the main nebular emission (v~158 km s _1 ) over a 
range Ad along the slit and divided from the data. A 
background level was then subtracted and the remaining 
counts within a wavelength (velocity) range A1-A2 were 
summed. The summed counts were converted into a to- 
tal flux. Surface brightness was determined by dividing 
the total flux by the area covered by the extraction region 
(Ad x slitwidth). Typical observed surface brightnesses 
were <~ 10~ 16 erg s _1 cm~ 2 arcsec~ 2 with a standard de- 
viation of 50%. 

There is a faint patch of Ha emission at 215 km s _1 
at the Position 1 LWRS aperture (Figure 3, top panel). 
This corresponds reasonably well with the observed O VI 
emission at 228 km s _1 . We measure the surface bright- 
ness of this emission patch to be 8.5 ± 2.4 x 10~ 17 
erg s _1 cm" 2 arcsec -2 over a slit range 57"-79" east of 
HD 5980. 



4. DISCUSSION 
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The detection of O VI emission as well as the kinemat- 
ics present in the Ha echelle data imply the presence of 
a strong shock in N66 associated with SNR 0057—7226. 
The SNR shock model was also favored by Hoopcs ct al. 
(2001) for the high velocity absorption towards HD 5980. 
Detection of the front side of SNR 0057—7226 has not been 
clearly established because the low velocity components 
cited by Koenigsberger et al. (2001) were detected at such 
a marginal level. In any event, the near side of the SNR 
is considerably more difficult to detect both in emission 
and absorption than the back side. Optical emission from 
the H II region complicates the situation and masks the 
morphology of the SNR. 

In the following discussion we adopt a working model of 
the physical association of various components in the N66 
region. A roughly spherical SNR is located on the near 
side of N66. The rear side of the SNR is propagating into 
relatively denser nebula while the near side is propagating 
through the more rarefied SMC ISM. HD 5980 lies behind 
the remnant, embedded within N66. The NGC 346 cluster 
stars lie outside the SNR, also embedded within N66. We 
determine the shock velocity, preshock hydrogen density 
and ram pressure of the shock wave using the emission 
and absorption sight lines toward SNR 0057— 7226. By ex- 
amining the surface brightness ratio between O VI and 
Ha, we measure the "completeness" of the shock. Then 
we discuss the relationship between the SNR and its sur- 
roundings and the implications of this analysis for observ- 
ing SNRs in other H II regions. 

4.1. Physical Properties 

In the N66 reference frame (vi sr =+158 km s _1 ), the 
material at Position 1 has a radial velocity of +70 km s _1 
while toward HD 5980, v r = +142 km s _1 . From X-ray 
and radio images, we choose the center of the remnant to 
be (J2000) a=00:59:27, £=-72:10:15, about 20" south of 
HD 5980. This position is also consistent with the location 
of the highest-velocity Ha emission (Figure 3). Using the 
projected distances from this center and the two radial ve- 
locity measurements, we derive a radius of curvature for 
the shock front of 80tio" an< ^ an expansion velocity for 
the O VI gas of v cxp = 147 ± 7 km s -1 . This places the 
observed X-ray edge on the eastern side of the shock front 
inclined ^60° to our line of sight. The expansion veloc- 
ity of 147 km s _1 is not the shock speed, but rather the 
bulk velocity of the postshock material which should be 
between 0.75 and 1.0 times the shock speed. 

First we caculate shock models for comparison against 
the observations. The models were calculated using an up- 
dated version of the code described by Raymond (1979). 
The main input parameters are the shock velocity, the 
preshock density and the elemental abundances. We ran 
models for shock velocities between 140 and 200 km s _1 
spaced by 10 km s _1 . A preshock hydrogen number den- 
sity of 1 cm~ 3 , and SMC elemental abundances (Russell & 
Dopita 1990) were used in all the models. The calculation 
is followed until the the recombination zone is complete 
and the gas temperature has reached about 1000 K. We 
present the intensities of selected lines from the shock mod- 
els in Table 3. The intensities are relative to Ha=100 and 
the flux from nearby multiplet lines of an ion have been 
summed. The O VI line flux produced by a shock increases 



sharply as a function of shock velocity in the range 150 - 
200 km s _1 . Shocks below ^160 km s" 1 do not produce 
much O VI and thus this represents a lower limit to the 
actual shock velocity in the SNR. Over this same range 
of shock velocities, the C III flux remains nearly constant. 
The changing ratio of 1(0 VI)/I(C III) is due to the in- 
creasing ionization of oxygen at higher velocities. 

As discussed in §3.1, the C III is strongly affected by 
absorption and scattering, and so it is not possible for 
us to derive an accurate assessment of the intrinsic ra- 
tio of 1(0 VI)/I(C III) in the spectrum of Position 1. 
However, we note that using the total (corrected) O VI 
flux with the observed C III flux places an upper limit of 
F(0 VI) /F(C III)<3. Ignoring the effect of differential red- 
dening for the moment, this clearly points us toward the 
lower end of the velocity range where O VI is produced. 
Based on this ratio, the observed expansion velocity and 
the shock velocity requirements for O VI production, we 
adopt 160 km s _1 as a representative shock velocity at 
Position 1 and use details from the model at this velocity 
below. 

We now compare the observed and calculated O VI 
surface brightnesses to estimate the preshock density re- 
quired to produce the observed O VI flux. The model 
predicts an O VI surface brightness (in both lines of the 
doublet) of I mo d = 6.08 x 10~ 5 erg s _1 cm~ 2 into 2ir stera- 
dians (for v s h O ck=160 km s _1 , n Q =l cm~ 3 , He/H=0.08) or 
I mod = 2.28 x 10~ 16 erg s _1 cm~ 2 arcscc~ 2 . In the opti- 
cally thin limit, intensity will scale with density so we can 
write 

lobs C1035 _ T 

l-mod Tlo 

where C1035 is the reddening correction factor at 1035 A, 
A is the surface area of shock (in arcseconds 2 ) within the 
aperture and n is the preshock density in units of cm~ 3 . 
lobs is the total O VI flux in both lines which we estimate 
by multiplying the corrected O VI A1032 flux by 1.5, thus 
lobs = 1-91 x 10~ 13 erg s _1 cm~ 2 . From Figure 1, we es- 
timate that ~| of the LWRS aperture is filled with X-ray 
emission. If the O VI emission follows the same pattern 
and the shock front is tilted 60° from our line of sight, 
then the aperture sampled emission from 600 arcseconds 2 
of shock front. 

Reddening is a potential source of uncertainty. Massey, 
Parker, & Garmany (1989) found the mean E(B— V)=0.14 
for the NGC 346 stellar population. However, since the 
SNR lies in front of the bulk of the cluster, this may 
overestimate the extinction to the SNR. On the other 
hand, using the spectral type and observed colors of 
Sk80, one can derive E(B-V)=0.11. We adopt this 
value and assume R v =3.1 (typical for the diffuse ISM) 
and the extinction curve of Hutchings & Giasson (2001) 
(E(1035A-V)/E(B-V)=18). This yields a reddening cor- 
rection Cio35=4.5. 

Solving the above equation, we get n = 6.3 cm" 3 at Po- 
sition 1, similar to an azimuthally- averaged electron den- 
sity of 7 cm~ 3 at the radius of Position 1 derived by Kcn- 
nicutt (1984). The dynamical pressure is then p v 2 s = 
1.3 m,H n a u 2 = 3.5 x 10 -9 dyne cm~ 2 . If more reddening 
is assumed and E(B— V) is raised to 0.14 (Cio35=6.8), the 
changes are minor: n Q = 9.5 cm~ 3 and P = 5.2 x 10 -9 
dyne cm~ 2 . Less reddening (E(B-V)=0.08, Ci 035 =3.0) 
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would yield n D = 4.2 cm 3 , P = 2.3 x 10 9 dyne cm 2 . 

This pressure is higher than those seen in regions of 
the Vela SNR. Sankrit et al. (2001) derive n G = 0.5 cm" 3 
and P=3.7xl0~ 10 dyne cm~ 2 for a X-ray bright, nearly 
face-on shock. Raymond et al. (1997) find P=1.6xl0 -9 
dyne cm~ 2 for an edge-on shock. Jenkins & Wallcrstcin 
(1995) find P=2-4xl0~ 10 dyne cm~ 2 for another locations 
in Vela while Kahn et al. '(1985) find P-lxlO" 10 dyne 
cm~ 2 for a region of diffuse X-ray emission. This vari- 
ation in pressures within the same remnant suggests lo- 
cal density and pressure enhancements caused by reverse 
shocks from blastwave-cloud interactions. We have only 
two widely-separated measurements for SNR 0057— 7226 
and can say nothing so specific but the high pressure and 
density at Position 1 as well as the X-ray emission mor- 
phology suggest similar shock-cloud interactions. 

The shock models discussed above predict the total 
column depths of ions produced in the post-shock flow. 
Hoopes et al. (2001) measured the column depths of sev- 
eral highly ionized species in the high-velocity component 
toward HD5980 (their Table 1) and conclude that a SNR 
shock is the most likely production mechanism for these 
ions in the observed quantities. Kocnigsberger et al. (2001) 
presented data on other ions not available to FUSE. In Fig- 
ure 5, we present the observed column depths of these FUV 
ions and the values predicted by the models with v s =160, 
180 and 200 km s _1 . The ions are ordered according to 
the ionization potential needed to produce them. The rel- 
ative column depths of all the ions down to Fe III are well 
matched by the 160 km s _1 model. The higher velocity 
models predict too much O VI relative to the other high 
ionization ions. The column depths scale linearly with 
prcshock density (all the models were run with no = 1 
cm -3 ). The offset between the 160 km s _1 model and 
the observations is 0.78 ± 0.14 dex, so a preshock density 
of n Q = 6 ± 2cm~ 3 can reproduce the observed column 
depths for the high ionization lines. The shock velocity 
and preshock density toward HD 5980 are thus in good 
agreement with our measurements for Position 1. 

The behavior of the low ionization lines is drastically dif- 
ferent. For Si II and Fe II, all the models predict column 
densities that arc an order of magnitude higher than the 
observations (see Figure 5). This dearth of low-ionization 
material can also be seen by comparing the observed and 
calculated ratio of O VI and Ha. The shock models (which 
have been calculated to the point of complete recombina- 
tion) predict 1(0 VI)/I(Ha)=3.3 for v s =160 km s" 1 and 
SMC abundances. The observed ratio at Position 1 is 
I{OV I) = SBjOVI) <7io35 /gq 
I(Ha) SB(Ha) C e563 fovi 

where SB indicates the surface brightness and / the aper- 
ture filling factor of the relevant emission. We assume the 
reddening correction Cio35=4.5 as discussed above; red- 
dening at Ha is moderate and we adopt C 6 563 = l-24. The 
filling factor / for each aperture is more uncertain. We 
assume fovi ~ § as discussed above. It is difficult to 
say where emission does and does not appear at a fine 
scale in the low S/N echelle data. The Ha aperture was 
chosen carefully to encompass only the faint emission in 
Figure 3, but we estimate that fua is still <~|. These 
values give 1(0 VI)/I(Ha)«18. Uncertainties in the fill- 
ing factors, reddening, and surface brightnesses probably 



make this ratio uncertain by a factor of 2, but the ratio is 
still significantly higher than the predicted ratio of 3.3. 

The weakness of Ha emission compared with O VI at 
Position 1 suggests the presence of an incomplete shock. 
There are many cases where incomplete shocks are ob- 
served in individual filaments or small portions of SNRs 
(Fesen, Blair, & Kirshncr 1982; Blair et al. 1991; Danforth, 
Blair, & Raymond 2001); however, these regions are all at 
sub-parsec scales. The LWRS aperture at Position 1 is 
sampling a region of SNR 0057— 7226 roughly 8 pc across. 
Furthermore, the observed Ha surface brightness does not 
vary strongly between different portions of the SNR im- 
plying that the conditions at Position 1 are typical of the 
remnant as a whole. In Figure 3, we see weak SNR Ha 
emission or none at all. This is a large object with mod- 
erate shock velocities implying that it is a middle-aged 
remnant like the Cygnus Loop or the Vela SNR; most of 
the optical emission from these remnants comes from com- 
plete, radiative shocks and these remnants are bright in 
Ha. If some portions of SNR 0057— 7226 were recombi- 
nationally complete, we would expect to see Ha emission 
<~10 times brighter than is observed. Thus, it appears that 
there is a global effect preventing the cooler parts of the 
recombination flow from forming. 

We argue that the gas at the back end of the post-shock 
flow is being prevented from recombining by the ionizing 
flux from the bright stars in N66. There are eleven stars 
of spectral type 06.5 or earlier located 2' (35 pc, in pro- 
jection) to the west of SNR 0057-7226. The combined 
ionizing flux from NGC 346 is 4.0xl0 40 erg s" 1 (Relano, 
Pcimbcrt & Bcckman 2002) or <~1.5xl0 51 ionizing pho- 
tons s' 1 . HD5980 itself contributes - 4 x 10 50 ionizing 
photons s _1 . Using basic ionization/rccombination equa- 
tions from Osterbrock (1989) we find that this photon flux 
in a medium with density ^5 cm~ 3 will yield an ionized 
sphere over 100 pc in radius. Indeed, Relano, Pcimbert & 
Bcckman (2002) find that N66 is density bounded and that 
45% of the ionizing photons from cluster stars escape to 
ionize the general SMC ISM. This zone of ionization could 
encompass SNR 0057— 7226 and affect the ionization bal- 
ance in the recombining post-shock gas. The dynamical 
time scale of the shock is of order a few hundred years. 
The photoionization time scales are much shorter (~100 
days), while the recombination times are longer (of order 
10,000 yr). Thus the ionizing flux will suppress Ha emis- 
sion. 

The strength and spectral energy distribution of the lo- 
cal ionizing flux will affect the column densities of the low 
and moderately ionized species. In theory, much of the 
discrepancy between the modeled and observed Fe II and 
Si II columns towards HD 5980 (Figure 5) could be due to 
iron and silicon being ionized to higher ionization stages. 
If the best fit model is scaled up to n Q = 6 cm~ 3 , the Fe III 
column is underpredicted by about 0.3 dex. This is clearly 
not enough to explain the 1.9 dex discrepancy in the Fe II 
column, but the stellar spectrum extends out to ~50 eV 
(Relano, Peimbcrt & Beckman 2002) and the local flux 
could be producing Fe IV. However, no such discrepancies 
are seen in the Si III and Si IV columns. 

A more likely scenario is that the sight line to HD 5980 
is passing through a small region where the swept up col- 
umn is lower than average for SNR 0057— 7226. While the 
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Fe III zone is complete, the Fe II zone would be incomplete 
regardless of ionizing radiation. A self-consistent shock 
model calculation that includes a photoionizing field is 
necessary to predict the column densities and the Ha flux 
accurately. While such modeling is beyond the scope of 
this paper, we conclude that the local ionizing field must 
have a significant effect on the observed characteristics of 
SNR 0057-7226 as a whole. 

4.2. The Relationship Between SNR 0057-7226 and N66 

The emission and absorption data discussed here lead 
to a schematic structure for the N66-SNR region shown 
in Figure 6. SNR 0057—7226 lies on the near side of N66 
and is encountering the denser nebular material on its back 
side. Position 1 is located at the edge of the X-ray emission 
(Figure 1). The observed radial velocity of O VI and Ha 
emission at Position 1 is ^225 km s _1 . This is lower than 
the observed velocity of absorption (+300 km s _1 ) toward 
HD 5980 which lies projected closer to the center of the 
remnant. Also it is ^75 km s _1 higher than the systemic 
velocity of N66 (+158 km s -1 ), which we would expect for 
material viewed tangentially at the edge of the expanding 
remnant. Either the systemic velocity of the SNR and N66 
differ by ~70 km s _1 , or the observed X-ray edge (a por- 
tion of which is Position 1) does not represent the actual 
limb of the SNR shock. The latter is assumed in Figure 6. 

We propose that the X-ray emission arises over a spheri- 
cal cap where the shock has encountered and heated denser 
material. The remnant is not yet even half-way submerged 
in the N66 material, and thus the observed radius of X- 
ray emission is smaller than the actual blast wave radius. 
The observed X-ray edge at Position 1 is inclined to our 
line of sight by -60°. HD 5980 and the other UV-bright 
stars are located within N66 and have ionized a region 
of surrounding material. The combined stellar winds and 
thermal pressure have formed a slowly expanding envelope 
around the cluster core. Since high-velocity absorption is 
not seen in the spectrum of Sk80 (Hoopes et al. 2001), it 
must lie within or in front of the SNR. HD 5980 lies be- 
hind the remnant and affects the ionization balance in the 
swept-up gas, preventing the formation of H I, S II, Fe II, 
Si II and other low-excitation species. The bright nebu- 
lar emission seen surrounding the UV bright stars of N66 
arises in dense material behind an ionized layer. 

The approaching SNR material is apparent in Ha in 
Figure 3, but has not been conclusively identified in FUV 
spectra. If the far side of the remnant is seen near ^300 
km s _1 and if the SNR is centered at the SMC systemic ve- 
locity of ~155 km s _1 , we expect the nearside absorption 
to be seen near v=10 km s . This region is essentially 
masked by Galactic absorption making it difficult to detect 
the SNR. Hoopes et al. (2001) compared the O VI spec- 
tra of Sk80 and HD 5980 near zero velocity and found no 
difference. However, both of these sight lines are toward 
the remnant and may contain nearly identical foreground 
O VI columns. 

The four FUSE targets in NGC 346 lie less than 2' from 
HD 5980 and are outside the SNR sight line. A compari- 
son of the averaged, normalized O VI profiles from these 
four targets and the two sight lines projected within the 
SNR is shown in Figure 7. We find a small but systematic 
excess in O VI absorption near v=0 km s _1 in the two 



remnant profiles. The difference represents an extra O VI 
column density of about 10 13 8 cm~ 2 or about j that of 
the backside value reported by Hoopes et al. (2001). We 
also find that the velocity centroid of the HD5980/Sk80 
zero velocity absorption is shifted by about +10 km s _1 
with respect to the NGC 346 Galactic absorption. This is 
not too different from the velocity found by Koenigsberger 
et al. (2001) for Si IV and C IV absorption features that 
they attribute to front-side SNR material. 

Recent work by Howk et al. (2002) on Galactic halo 
O VI absorption has shown significant variations in O VI 
column over as little as 3'. The variation in N(0 VI) be- 
tween the four NGC 346 targets (mean separation 19") is 
9% (the maximum variation between two measurements 
is 12%) while the variation between the two SNR sight 
lines (separation <1') is only 8%, comparable to the un- 
certainties due to continuum placement, noise, and limits 
of integration. However, the cluster and SNR stars show a 
relative difference of 53% in the zero- velocity O VI column 
over ~2'. Because of the potential Galactic halo variation 
on small scales, the matter of frontside absorption remains 
inconclusive. However, it is plausible we are seeing front 
side absorption near v=+10 km s _1 in the HD 5980/Sk80 
data. If true, then Sk80 would be inside SNR 0057-7226, 
as shown in Figure 6. Even if it is not, we can place an 
upper limit of N(0 VI) < 6 x 10 13 cm~ 2 for the frontside 
O VI column density, much less than is observed from the 
back side of SNR 0057-7226. 

In a recent paper, Velazquez, Koenigsberger, & Raga 
(2003) model the collision of a SNR shock with a wind- 
blown bubble around an evolved star and suggest that the 
far side of SNR 0057— 7226 is interacting with the stellar 
wind of HD 5980. The X-ray image of the system (Naze 
et al. 2002) does, in fact, show some resemblance to the 
models and there is reasonable agreement between pre- 
dicted and observed X-ray fluxes. Our observations allow 
us to further test this wind-SNR model. The Velazquez 
et al. wind-blown bubble has a radius of ~10 pc while our 
Position 1 (at the rim of the SNR) is at least twice that 
distance (in projection) from HD 5980; thus, Position 1 is 
sampling the ISM, not stellar wind. In §4.1 we derived a 
preshock density of ~6 cm' 3 for the material at both Po- 
sition 1 and in the shock toward HD 5980, at least two or- 
ders of magnitude higher than the interior of the modeled 
wind bubbles (see Figures 1 and 2 of Velazquez, Koenigs- 
berger, & Raga 2003). Furthermore, the +300 km s _1 
component exhibits an ionization structure expected for a 
single SNR shock (Hoopes et al. 2001). If this component 
arose in a turbulent region between SNR and stellar wind 
shock fronts, as suggested by Velazquez, Koenigsberger, & 
Raga (2003), we might not expect to see the kinematics 
structure observed by Hoopes et al. (2001). Clearly, higher 
resolution spectral observations and detailed shock models 
are required. 

4.3. Implications for SNRs in H II Regions 

Core-collapse supernovae arise from massive stars often 
formed in clusters and associations-the same massive stars 
which photoionize their surrounding gas and create bright 
H II regions. Massive SN progenitors cannot migrate far 
from their birth places in their short lifetimes and thus 
SNRs associated with H II regions should occur frequently. 
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To various degrees, this is not seen in optical SNR surveys 
of nearby galaxies. 

The primary optical method for detecting SNRs, es- 
pecially in galaxies beyond the Magellanic Clouds, is to 
compare Ha to [S II] emission, and pick out objects with 
[S Il]/Ha> 0.4 (cf. Long ct al. 1990; Blair & Long 1997, 
and references therein). In recent years, as the sensitivity 
and resolution of radio and X-ray observations has become 
competitive with optical ground-based observations, inde- 
pendent searches at these wavelengths have been under- 
taken and compared against the optical SNR lists (Lacey, 
Duric, & Goss 1997; Gordon ct al. 1998, 1999; Pannuti 
et al. 2000, 2002). These multiwavelcngth surveys and 
comparisons provide a more complete survey method, but 
in general do not find many of the same SNRs. Discus- 
sions of selection effects usually point toward "confusion" 
as being responsible for incompleteness of optical SNR sur- 
veys in crowded and complex regions of emission. Clearly 
this is plausible at some level, but the denser conditions in 
H II regions might also be expected to cause brighter op- 
tical (radiative) SNRs since the shock emissivity increases 
as the density squared. Scanning through the optical SNR 
catalogs of nearby galaxies (Blair & Long 1997; Long et 
al. 1990; Matonick & Fesen 1997, to name a few), one can 
find a number of stunning counter examples to the "con- 
fusion" hypothesis. There are often bright optical SNRs 
found in confused regions and adjacent to bright H II re- 
gion emission. The mechanism we propose in §4.2 seems to 
be a more plausible explanation for why the optical SNR 
searches miss some SNRs and not others in confused re- 
gions. If a supernova occurs close enough to ionizing stars, 
the potential bright optical SNR emission in Ha and S II 
is prevented from forming and no optical SNR can be de- 
tected, especially against bright nebular emission. 

Our investigation of SNR 0057— 7226, however, demon- 
strates that bright FUV emission lines are still present 
in such a situation and may provide an alternate way of 
not only detecting but studying such SNRs in detail. De- 
tecting SNRs with higher ions circumvents the problems 
of truncated recombination and bright background emis- 
sion. Stellar photoionization will not prevent the produc- 
tion of high ions nor will they be contaminated signifi- 
cantly by surrounding nebular emission. Our observations 
have shown that in a bright H II region, the detection of lo- 
calized emission from high-ionization species such as O VI 
can provide a clear signature of SNR emission. Position 1 
lies on the edge of the SNR and shows strong emission in 
O VI and C III. Position 2 lies outside the remnant but 
still within bright optical emission; yet no line emission is 
present. In conjunction with radio and X-ray data, such 
a detection is conclusive. We note however that the use 
of the FUV is restricted to lines of sight where extinction 
is low. Nearby extragalactic targets such as the Magel- 
lanic Clouds and Local Group galaxies will be the easiest 



to investigate with this technique. 

5. SUMMARY 

We see strong emission in O VI and C III at a position 
on the eastern X-ray edge of SNR 0057-7226. Absorp- 
tion in these same ions is seen at ^300 km s^ 1 toward the 
bright LBV system HD 5980 which lies behind the rem- 
nant. Another position away from the remnant but within 
the bright nebular extent of N66 shows no FUV line emis- 
sion. The remnant lies on the near side of N66 and the 
O VI is formed in shocks as the blast wave encounters the 
denser material of the nebula. 

Seven longslit echelle spectra in Ha and [N II] show 
that N66 has turbulence similar to other Magellanic H II 
regions. Two slowly expanding features are seen: the com- 
pact H II region N66C and a bright shell around the core of 
NGC346. Three of the echelle spectra show faint patches 
of Ha emission at high and low velocities which appear 
as a rough expansion feature corresponding spatially and 
kinematically to SNR 0057-7226. 

Observed C III and O VI flux compared with shock mod- 
els allows us to estimate the shock velocity and determine 
the ambient density and ram pressure at Position 1. We 
use observed column depths from Hoopes et al. (2001) and 
Kocnigsberger ct al. (2001) to derive the shock velocity 
and density of the SNR shock toward HD 5980. These val- 
ues are in good agreement with each other and with other 
values in the literature. 

Weak Ha emission as well as low column densities of 
low-ionization species such as Fe II show that the ioniza- 
tion balance in the post-shock gas is affected by ionizing 
stellar photons. The ionizing flux from the many hot stars 
in NGC 346 is sufficient to keep swept-up material out of 
ionization equilibriuim. This highlights the importance 
and inherent difficulty of looking for SNRs in bright H II 
regions. Direct optical identification is complicated by the 
bright nebular emission. FUV emission and absorption, 
along with X-ray and radio observations, are good ways to 
circumvent these difficulties. 
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Table 1 
FUSE Observation Summary 



Sight Line 


RA (J2000) 
h m s 


Dec (J2000) 
/ // 


SpType 


FUSEID 


Aper a 


#Exp 


ExpTime 
(sec) 


Obs Date 


NGC 346- WB6 


00 58 57.74 


-72 10 33.6 


04V((f)) 


P20305 b 


MDRS 


5 


10992 


2001-Sep-25 


NGC 346- WB4 


00 59 00.39 


-72 10 37.9 


05-6V 


P20304 b 


MDRS 


6 


11853 


2001-Aug-25 


NGC 346- WB3 


00 59 01.09 


-72 10 28.2 


02III(f*) 


P20303 b 


MDRS 


3 


8482 


2001-Aug-25 


NGC 346- WB1 


00 59 04.81 


-72 10 24.8 


04III(n)(f) 


P20302 b 


MDRS 


3 


4625 


2001-Aug-25 


HD 5980 


00 59 26.55 


-72 09 53.8 


WN var 


P10301 


LWRS 


4 


5734 


2000-Jul-02 


Sk80 


00 59 31.95 


-72 10 45.8 


07 Iaf+ 


P10302 


LWRS 


4 


11699 


2000-Jul-02 


N66-Pos 1 


00 59 42 


-72 09 49 




P20305 


LWRS 


5 


10992 


2001-Sep-25 


N66-Pos 2 


00 59 37 


-72 07 57 




P20302 


LWRS 


3 


4625 


2001-Aug-25 



a FUSE Apertures: LWRS=30" square, MDRS=4"x20". 
b Combined into composite NGC 346 spectrum. 
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Table 2 

Summary of CTIO 4-m Observations 



Target RA (J2000) Dec (J2000) PA ExpTime Date 

h m s o , „ o sec y m d 



HD 5980 NS 


00 


59 


26.55 


-72 


09 


53.8 





1200 


2001-Dec-24 


HD 5980 EW 


00 


59 


26.55 


-72 


09 


53.8 


90 


1200 


2001-Dec-25 


Sk80 NS 


00 


59 


31.95 


-72 


10 


45.8 





1200 


2001-Dec-26 


Sk80 EW 


00 


59 


31.95 


-72 


10 


45.8 


90 


600 


2000-Dcc-lO 


NGC 346- WB6 


00 


58 


57.74 


-72 


10 


33.6 


90 


600 


2000-Dec-10 


NGC 346-WB3+4 


00 


59 


01.09 


-72 


10 


28.2 


197 


600 


2000-Dec-10 


NGC 346- WB1 


00 


59 


04.81 


-72 


10 


24.8 





600 


2001-Dec-26 



Table 3 

Selected Lines from the Shock Models 



Line / v s (km s 1 ) 


140 


150 


160 


170 


180 


190 


200 


Ha 6563 


100 


100 


100 


100 


100 


100 


100 


S VI 937 


12 


19 


18 


20 


23 


27 


30 


C III 977 


352 


304 


275 


258 


247 


239 


232 


N III 991 


17 


13 


12 


11 


10 


10 


10 


Ne VI 1006 





1 


4 


12 


27 


45 


58 


O VI 1034 


9 


90 


334 


708 


996 


1156 


1236 


S IV 1070 


11 


8 


8 


7 


7 


7 


7 


Ne V 1146 


3 


8 


13 


18 


24 


27 


28 


N V 1240 


17 


27 


23 


23 


24 


24 


24 


Si IV 1397 


31 


28 


26 


25 


24 


24 


23 


O IV] 1402 


100 


112 


108 


101 


97 


94 


92 


N IV] 1490 


9 


7 


6 


5 


5 


5 


5 


C IV 1549 


211 


174 


156 


148 


145 


142 


140 


O III] 1662 


59 


52 


47 


44 


42 


41 


40 


N III] 1748 


6 


5 


4 


4 


4 


4 


4 


I(Ha) a 


1.268 


1.551 


1.819 


2.048 


2.249 


2.431 


2.622 



a Flux units are 10 5 erg s 1 cm 2 into 2tt steradian for n =l cm 3 and SMC 
abundances. 
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Fig. 1. — The giant H II region N66. An Ha image is shown in the left panel and FUSE stellar targets are marked with small squares. 
LWRS emission apertures are the large boxes in both panels. Dashed lines show the locations of the echelle spectral slits. The lowest X-ray 
contour is also shown. The right panel shows the Chandra 0.25-10 keV X-rays (courtesy of Naze et al. 2002). The ACIS chip boundaries can 
be seen partially obscuring the southern edge of the remnant. 



11 



T 3 

•< 

'to 2 

e i 



C -1 



cm 



Posl 



-2 
965 



Pos2 



970 



975 980 
Wavelength (A) 



985 



m 2 

CM 

o 1 



FT "2 



Posl 



Pos2 




1020 



1025 



1030 1035 
Wavelength (A) 



1040 



1045 



Fig. 2.— FUSE emission spectra at Positions 1 (SNR 0057-7226) and 2 (N66 background) around 977 and 1032A. Position 2 data has 
been offset downward by two units. Note the strong O VI and C III emission from Position 1 and the lack of emission from Position 2. 
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Fig. 3.— Ha 6562.82A echelle observations of SNR 0057-7226. The top two slits were centered on HD 5980 oriented cast-west (top) and 
north-south (middle). The bottom slit shows a north-south slit centered on Sk80. Bright nebular emission is seen at v~ +160 km s _1 and 
the stellar continuum from the indicated star is a vertical band in the center. Faint, patchy emission at high and low velocities can be seen 
above and below the systemic velocity along all three slits; this emission arises in the SNR shock. Terrestrial Ha airglow can be seen as a 
faint horizontal stripe at v~ —20 km s _1 . The location of the Position 1 LWRS slit is shown by the dashed lines in the top panel. 
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Fig. 4. — The VI A1032 emission (top panel) is essentially free from H2 absorption and is affected only on the blue side by foreground 
Galactic and SMC O VI absorption. The emission has been fit as described in the text (solid curve). The intrinsic Gaussian implied by the fit 
is shown as the dashed curve. We scale this intrinsic profile by a factor of two in the middle panel (dashed curve) and quantify the absorption 
in the more heavily contaminated O VI A1038 line. The solid line represents the model fit. The fits are consistent with a 2:1 optically thin 
ratio for the O VI lines. The bottom panel shows C III emission along with the O VI A1032 intrinsic profile for comparison. No fitting has 
been attempted on the C III data. 
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Fig. 5. — Observed and modeled absorption columns for various ions toward HD5980. The squares and diamonds show measured column 
density and quoted uncertainties from Hoopes et al. (2001) (squares) and Koenigsberger et al. (2001) (diamonds). Ions are shown in order 
of decreasing ionization potential from O VI to Fe II. Measured columns for S VI, Fe III, and Fe II are upper limits. The dotted lines show 
the columns predicted by shock models. In all three models, SMC abundance and preshock density n = 1 cm" 3 are used. The observed 
columns are offset from the 160 km s -1 model by 0.78 ± 0.14 dex. Since the model columns scale linearly with preshock density, this gives 
n Q = 6 ± 2 cm~ 3 . 




Fig. 6. — A schematic representation of the proposed physical relationship between SNR 0057— 7226 and N66 is shown. Our line of sight 
is assumed to be from the top of the Figure. SNR 0057— 7226 is located on the near side of N66 and is interacting with the denser, ionized 
gas on the rear side. The O VI, C III, and X-ray emission seen arises from this shock interaction. The approaching side of the shock is 
propagating through a more diffuse medium and is harder to detect. HD 5980 is located behind the SNR while Sk80 is either within or in 
front of the SNR. The massive stars in the field have ionized a layer of nebular material and keep the swept-up gas behind the Shockwave 
from recombining completely. A slowly expanding bubble around NGC 346 is powered by stellar winds and thermal pressure. 
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Fig. 7. — The normalized O VI absorption from HD 5980 (dashed) and Sk80 (dotted) are markedly different from the mean NGC346 data 
(solid) near v=0 km s . Wc propose that this excess absorption represents the approaching (blue-shifted) side of the SNR. Absorption from 
the receding side is clearly seen in the spectrum of HD5980 near 300 km s _1 , as reported by Hoopcs et al. (2001) but is not evident in the 
other spectra. 



